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Abstract 

Conventional thermo-statistics address infinite homogeneous systems within the 
canonical ensemble. (Only in this case this is equivalent to the fundamental mi- 
crocanonical ensemble.) However, some 170 years ago the original motivation of 
thermodynamics was the description of steam engines, i.e. boiling water. Its essen- 
tial physics is the separation of the gas phase from the liquid. Of course, boiling water 
is inhomogeneous and as such cannot be treated by conventional thermo-statistics. 
Then it is not astonishing, that a phase transition of first order is signaled canoni- 
cally by a Yang-Lee singularity. Thus it is only treated correctly by microcanonical 
Boltzmann- Planck statistics. It turns out that the Boltzmann-Planck statistics is 
much richer and gives fundamental insight into statistical mechanics and especially 
into entropy. This can be done to a far extend rigorously and analytically. As no 
extensivity, no thermodynamic limit, no concavity, no homogeneity is needed, it 
also applies to astro-physical systems. The deep and essential difference between 
"extensive" and "intensive" control parameters, i.e. microcanonical and canonical 
statistics, is exemplified by rotating, self-gravitating systems. In the present paper 
the necessary appearance of a convex entropy S{E) and negative heat capacity at 
phase separation in small as well macroscopic systems independently of the range 
of the force is pointed out. Thus the old puzzle of stellar statistics is finally solved, 
the appearance of negative heat capacity which is forbidden and cannot appear in 
the canonical formalism. 
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1 Introduction 

Conventional statistical mechanics addresses homogeneous macroscopic sys- 
tems in the thermodynamic limit. These are traditionally treated in canonical 

ensembles controlled by intensive temperature T, chemical potential fi and/or 
pressure P. In the canonical ensemble the heat capacity is given by the fluc- 
tuation of the energy <{6E)'^> /T^ > 0. 

As in astro-physics the heat capacity is often negative it is immediately clear 
that astro-physical systems are not in the canonical ensemble. This was often 
considered as a paradoxical feature of the statistics of self-gravitating sys- 
tems. Here we will show that this is not a mistake of equilibrium statistics 
when applied to self-gravitating systems but is a generic feature of statistical 
mechanics of any many-body systems at phase separation, independently of 
the range of the interactions, ref.(l). 

As the original motivation of thermodynamics was the understanding of boil- 
ing water in steam-engines, this points to a basic misconception of conventional 
canonical thermo-statistics. As additional benefit of our reformulation of the 
basics of statistical mechanics by microcanonical statistics there is a rather 
simple interpretation of entropy, the characteristic entity of thermodynamics. 



2 What is entropy? 

Boltzmann, ref.(2), defined the entropy of an isolated system in terms of the 
sum of all possible configurations, W, which the system can assume consistent 
with its constraints of given energy, volume, and further conserved constraints: 



in semi-classical approximation. E is the total energy, is the number of par- 
ticles and V the volume. Or, more appropriate for a finite quantum-mechanical 



S=k*lnW 




as written on Boltzmann's tomb-stone, with 





(D.H.E.Gross). 
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system: 



W{ E, TV, V) = Tr[VE] = E ^ ^ "^^^^ ^^"^''^ (3) 

and < K < + eo 



and eo ~ the macroscopic energy resolution. This is still up to day the deep- 
est, most fundamental, and most simple definition of entropy. There is no 
need of the thermodynamic limit, no need of concavity, extensivity, and homo- 
geneity. Schrodinger was wrong saying that microcanonical statistics is only 
good for diluted systems, ref.(3). It may very well also address the sohd-liquid 
transition ref.(4) and even self- gravitating systems as we will demonstrate in 
this article. In its semi-classical approximation, eq.(2), W{E, N,V,-- ■) simply 
measures the area of the sub-manifold of points in the 6A^-dimensional phase- 
space (F-space) with prescribed energy E, particle number N, volume V, and 
some other time invariant constraints which are here suppressed for simplicity. 
Because it was Planck who coined it in this mathematical form, I will call it 
the Boltzmann-Planck principle. 

The Boltzmann-Planck formula has a simple but deep physical interpretation: 
1^ or 5" measure our ignorance about the complete set of initial values for all 
6A^ microscopic degrees of freedom which are needed to specify the A^-body 
system unambiguously, ref.(5). To have complete knowledge of the system we 
would need to know [within its semiclassical approximation (2)] the initial 
positions and velocities of all N particles in the system, which means we 
would need to know a total of 6A^ values. Then W would be equal to one 
and the entropy, S, would be zero. However, we usually only know the value 
of a few parameters that are conserved or change slowly with time, such as 
the energy, number of particles, volume and so on. We generally know very 
httle about the positions and velocities of the particles. The manifold of all 
these points in the 6A^-dim. phase space, consistent with the given conserved 
macroscopic constraints of E,N,V,- ■ ■, is the microcanonical ensemble, which 
has a well-defined geometrical size W and, by equation (1), a non- vanishing 
entropy, S{E, N,V,---). The dependence of S{E, N,V,-- •) on its arguments 
determines completely thermostatics and equilibrium thermodynamics. 

Clearly, Hamiltonian (Liouvillean) dynamics of the system cannot create the 
missing information about the initial values - i.e. the entropy S{E, N,V,-- ■) 
cannot decrease. As has been further worked out in ref.(6) and more recently in 
ref.(7) the inherent finite resolution of the macroscopic description implies an 
increase ofW or S with time when an external constraint is relaxed, c.f. chapter 
5. Such is a statement of the second law of thermodynamics, ref.(8), which 
requires that the internal production of entropy be positive or zero for every 
spontaneous process. Analysis of the consequences of the second law by the 
microcanonical ensemble is appropriate because, in an isolated system (which 
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is the one relevant for the microcanonical ensemble), the changes in total 
entropy must represent the internal production of entropy, sec above, and 
there are no additional uncontrolled fluctuating energy exchanges with the 
environment. 



3 No phase separation (no boiling water) without convex non- 
extensive S{E). 

The weight e'^^^)"^/-^ of configurations with energy E in the definition of the 
canonical partition sum 



becomes here bimodal, at the transition temperature it has two peaks, the 
liquid and the gas configurations which are separated in energy by the latent 
heat. Consequently S{E) must be convex (hke y — x'^) and the weight in (4) 
has a minimum between the two pure phases. Of course, the minimum can 
only be seen in the microcanonical ensemble where the energy is controlled and 
its fluctuations forbidden. Otherwise, the system would fluctuate between the 
two pure phases (inter-phase fluctuation) by an, for macroscopic systems even 
macroscopic, energy AE ~ Eiat <x N oi the order of the latent heat. I.e. the 
convexity of S{E) is the generic signal of a phase transition of first order and of 
phase- separation, ref.(9). This applies also to macroscopic systems coupled by 
short or long-range interactions. Of course, even there phase-separation exists 
and there is a minimum (convexity) in S{E) and, consequently, negative heat 
capacity(l). 

It must be emphasized: In conventional canonical statistics phase-separated 
systems are treated by artificial constructions of two coexisting simple configu- 
rations of each phase hke a single drop with vapor around (Wulff-construction) . 
But this is not a common ensemble of various droplets with its inter-phase 
fluctuations leading e.g. to a negative heat capacity. Such macroscopic en- 
ergy fluctuations and the resulting negative heat capacity are already early 
discussed in high-energy physics by Carlitz, ref.(lO). 



4 AppUcation to Eistrophysics 

Classical Boltzmann-Planck microcanonical statistics applies very well to the 
equilibrium of self-gravitating systems provided we put the system into a box 



oo 




(4) 
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Fig. 1. Convex intruder in the entropy as function of an conserved order parameter. 
Gibbs double tangent rule = concave hull. The derivative gives the Maxwell line 

in T{E). Three values of the order parameter correspond to the single intensive 
"temperature" (inverse slope T = 1/A). At oi the system is in the "liquid" phase, 
at 03 in its "gas" phase. The difference 03— oi is the latent "heat" and of macroscopic 
size oc N. 

(ignoring slow evaporation) and use a hard core in the interaction at short 
distances where anyhow non-gravitational physics like hydrogen burning dom- 
inates gravity. Padmanabhan in his stimulating review ref.(ll) discussed this 
point in detail. 

We simplify the original many-body problem by the mean-field approxima- 
tion, details in ref(7). Here the (certainly also important) many-body corre- 
lations are ignored. The entropy is then a function of the one-body densities 
p(r). S{E,p) is maximized and a closed equation for the optimal one-body 
density distributions obtained. To get this result and to get finite p(r) the 
short-distance cut-off and a hard core in the two-body interaction U{t-t') is 
essential. The resulting mean- field theory therefore does not lead to an isother- 
mal sphere. The latter is a solution of Poisson's equation (ref.(ll)) with a strict 
1/r interaction. 

Moreover, for equilibrium one does not need any exotic q-deformed statis- 
tics, ref.(13; 14). The error of rcf.(14) is the use of a generalized variant of 
Boltzmann-Gibbs canonical ensemble instead of the basic microcanonical one, 
a point on which Tsallis seems to agree in his discussion in ref.(14), see also 



In fig. 2 we show the equilibrium distributions of self-gravitating systems under 
large angular-momentum in the mean-field approximation: 

The necessity of using "extensive" instead of "intensive" control parameters 



ref.(7). 
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Fig. 2. Rotating multi-star-systems at equilibrium. The left shows a rotating dou- 
ble-star system in the DC phase, c.f. fig. 3. This is an inhomogeneous phase, analog 
to nuclear multifragmentation. The right is one of the many different, unstable, con- 
figurations existing in the mixed phase with negative heat capacity. Here the system 
fluctuates between such ring systems, systems of stars rotating around a central star 
but also mono-stars and eventual gas. This region is very interesting but must still 
be more investigated. [Cover page of Phys.Rev.Lett. vol 89, (July 2002)] 

is explicit in astrophysical problems. E.g.: for the description of rotating stars 
one conventionally works at a given temperature and fixed angular velocity 
Q c.f. ref.(15). Of course in reality there is neither a heat bath nor a rotating 
disk. Moreover, the latter scenario is fundamentally wrong as at the periphery 
of the disk the rotational velocity may even become larger than velocity of 
light. Non-extensive systems like astro-physical ones do not allow a "field- 
theoretical" description controlled by intensive fields ! 

E.g. configurations with a maximum of random energy on a rotating disk, i.e. 
at fixed rotational velocity Q (which is intensive): 



random 



E 



-E. 



pot 



(5) 



and consequently with the largest entropy are the ones with smallest moment 
of inertia O, compact single stars. Just the opposite happens when the angular- 
momentum L (extensive) and not the angular velocity Q is fixed: 

^2 

Erandom E „ EpQf. (6) 



Then configurations with large moment of inertia are maximizing the phase 
space and the entropy. I.e. eventually double or multi stars are produced, as 
observed in reality. 

In figure 3 one clearly sees the rich and realistic microcanonical phase-diagram 
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Fig. 3. Phase diagram of rotating self-gravitating systems in the en- 
ergy-angular-momentum L)-plane. DC: region of double-stars, G: gas phase, 
SC: single stars. In the mixed region one finds various exotic configurations like 
ring-systems in coexistence with gas, double stars or single stars. In this region 
of phase-separation the heat capacity is negative and the entropy is convex. The 
dashed lines E — L = —1 (left) and E = L (right) delimit the region where sys- 
tematic calculations were carried out. At a few points outside of this strip some 
calculations like the left of fig. (2) were also done. 

of a rotating gravitating system controlled by the conserved "extensive" pa- 
rameters energy and angular- momentum, ref.(16). 



5 On the Second Law 



Also the Second Law becomes significantly more transparent within micro- 
canonical statistics: 

Many formulations of the second law exist. The understanding of entropy is 
sometimes obscured by frequent use of the Boltzmann-Gibbs canonical ensem- 
ble, and thermodynamic limit; and the relationship of entropy to the second 
law is often beset with confusion between external transfers of entropy and its 
internal production. Perhaps the clearest statement is by Prigogine, ref.(8): 
He distinguishes between internal and external entropy transfer. Whereas the 
transfer df.,j.tS = dQ/T to an external heat bath can be positive, negative, or 
zero, the internal generation of entropy dintS within the system is necessarily 
dintS > 0. This is the sharpest formulation of the second law. dintS is best 
seen in the microcanonical ensemble, as there are no couplings to any heat 
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bath and no sometimes uncontrolled exchanges of energy with it. 



The second law is then most rigorously given by equation (7): 



dS — dextS + dintS 
dextS = positive, negative, or zero 

d^ntS>0 (7) 



A geometric visualization how an equilibrium distribution in the 6iV-dim phase 
space expands into a (fat)-fractal (Gibbs ink lines) after a constraint, here the 
wall between the phase-space volumes between Vi and V2 is suddenly released 
is given by figure (4). 




^ < ^0 t > to 

Fig. 4. The sudden expansion in phase-space of an equilibrized system left into a 
larger phase-space at t > to right. When at to the wall closing volume Va is suddenly 
removed, the compact set M{t < Iq), left side, develops into an increasingly folded 
"spaghetti" -like, though still compact, distribution in phase-space with rising time 
t > to, the Gibbs' "ink-lines". The right figure shows only the early form of the 
distribution. At much larger times it will become more and more fractal and finally 
dense in the new phase space. The grid illustrates the boxes of the box-counting 
method. All boxes with size 6^^ which overlap with M(t) are counted by Ns in 
eq.(8) 



The Liouville theorem tells us that the (Riemannian-) size of it as given by 
eq.(2) does not change. Keeping in mind that due to the redundant information 
about the system (c.f. chapter 2) it is in general impossible to distinguish a 
point inside the fractal from a neighboring point outside. Therefore the proper 
way to "measure" the size of the manifold is by its box-counting "measure" . 
I.e. the phase-space integrals in eq.(2) or eq.(3) have to be calculated by the 
box-counting volume volbox^ as is explained in fig. (4) (details in refs.(6; 7)): 

yohox[M{E, N, t > to)] := Mls^o^''' Ns[M{E , N, t > to)] (8) 
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which is > yolbox[^{E, N,t < to)]. That is the Second Law. 



6 Further challenges 

The story of applying microcanonical statistics to self-gravitating systems and 
eventually to cosmology is far open. There are many questions to answer, e.g.: 

There is the problem of statistics in an expanding universe as discussed by 
Padmanabhan ref.(ll), the question of defining Boltzmann-Planck statistics 
in general relativity and connected the problem of avoiding field-theoretical 
descriptions with intensive statistical control-parameters or using them despite 
the fact that first order transitions demand extensive control parameters as 
discussed here. The present author, coming from far away from nuclear and 
statistical physics, is open for any helpful suggestion. An interdisciplinary 
meeting like this NBS2005 is here very important. 1 like to thank Prancoise 
Combes and Raoul Robert for this great opportunity. 



7 Conclusion 

Conventional canonical statistics is insufficient to handle the original goal 
of Thermodynamics, phase separations with their characteristic inter-phase 
fluctuations. Here the microcanonical ensemble is much richer and leads to 
negative specific heat. This is thus not a characteristics of self-gravitating sys- 
tems. Phase transitions can be found without the thermodynamic limit, without 
homogeneity, concavity, and extensivity. Microcanonical statistics opens to a 
much wider class of applications as e.g. astrophysical systems. Moreover it 
gives a simple formulation and proof of the second law. 
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